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Abstract 

A ID model hamiltonian that is motivated by the recent discovery of the heavy- fermion 

behavior in the cuprates of the Nd2Cu04 type is studied. It consists oit — J interacting 

conduction electrons coupled to a lattice of localized spins through a Kondo exchange 

term Jk- Exact diagonalization and density matrix renormalization group methods are 

used. The latter method is generalized to arbitrary densities. At half-filling, a spin gap 

opens for all Jk > 0. Away from half-filling, it is shown that, at strong Jk, the ground 

state is an unsaturated ferromagnet . At weak Jk the system is in a paramagnetic 

phase with enhanced RKKY correlations. The importance of self-screening of the local 

moments in the depletion regime is discussed. We argue that these findings transcend 

the specifics of the model. 

PACS: 75.30.Mb, 75.40.Mg, 71.27.+a,75.20.Hr 
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The study of quantum moments embedded in a metallic host is one of the most active field 
of condensed matter physics. The main problem is the understanding of the low-temperature 
behavior of localized moments, which are typically rare-earth lanthanide or actinide ions, 
interacting with p or d bands of conduction electrons. It has been argued that there are 
two competing effects: the quenching of the local spins through the Kondo effect and their 
magnetic ordering due to the RKKY exchange interaction The canonical Hamiltonians 
of heavy-fermion compounds are the periodic Anderson lattice model (PAM) and the Kondo 
lattice model (KLM). The PAM and the KLM describe a lattice of local moments coupled 
to a wide conduction band. It is widely accepted that this situation prevails in the so-called 
heavy-fermion materials 0. 

The recent discovery of a heavy-fermion phenomenon in the Ce-doped Neodynium 
cuprate Ndi^CcQ ^CuO^ has led to an increasing interest in the study of strongly corre- 
lated electrons coupled to magnetic moments [Q. Although these compounds are at least 
two-dimensional, as a first step towards their understanding, we report a study of a ID 
Hamiltonian which retains the basic ingredients of their physics. This Hamiltonian de- 
scribes the interplay between strongly correlated electrons and a lattice of local moments. 
The strongly correlated conduction electron system is described by the ID t — J model. The 
interaction between the conduction electrons and the Nd ions is represented by a Kondo 
exchange term. Using standard notations, the Hamiltonian is written as follows: 

iu 

+ -^Jf Xl(^jc ■ Sj+ic ^^J^) + -^i^ XI ^jc ■ Sj/ (1) 

i i 

The Hamiltonian (1) is a natural extension of the one- dimensional KLM in the limit 
of strongly correlated electrons. We recall that double occupancy is forbidden here. The 
indices c and / stand for conduction and localized electrons respectively. We set the hopping 
parameter t and the Heisenberg coupling Ju equal to 1. For those parameters, which are 
relevant for the cuprates, the ground state of the \D t — J model is a spin density wave . 
We discuss the nature of the ground state as the band filling p of the conduction electrons 

2 



and Jk are varied. 

We use the density matrix renormalization group {DMRG) method of White 0. The 
calculations are done using the infinite system method with 3 blocks instead of four as usual. 
One can immediately realize that there is a problem to fix the electron density at a given 
value at each iteration. To overcome this, we have targeted the two states whose electron 
number brackets the desired density when building the density matrix, This algorithm has 
been successfully checked for the one- dimensional t — J model 0. In the present study, we 
have kept up to m = 140 states in the two external blocks. The states were labelled by the 
electron number and by the total 2;-component of the spin. The maximum lattice size we 
have reached is 75 . The truncation error defined as 1 —p{m), where p{m) is the sum of the 
eigenvalues of the density matrix of the states kept, was of the order of 5 x 10~^. 

At half-filling, the t — J chain reduces to a Heisenberg antiferromagnet. We have found 
that the singlet ground state is separated from the first excited triplet state by a gap for 
all Jk > 0. As shown in Fig.l, the Jk dependence of this gap A is very close to that 
of the Kondo necklace model (KNM) It is linear in the strong Jk limit, and, has 
probably a non analytic form for small Jk- We note that the values of the gap found here 
in the weak coupling limit are much larger than that in the KNM. Thus, the interaction 
between the conduction electrons, represented here by the z part of the Heisenberg term 
after the Jordan- Wigner transformations, tends to enhance the gap. In order to understand 
the character of the gap, we have compared it with the binding energy per site which is 
defined as Eb = [Eg{Jk = 0, A^) — Eg{Jk, N)] /N, where Eg{Jk,N) is the ground state 
energy. It can be observed in Fig.l that in the strong coupling limit, A is larger than 
Eb'- the excitations of the system are obtained by breaking up on-site singlets. In the 
weak coupling regime, Eb becomes greater than A, this means that the system can be 
excited without the destruction of local singlets. The gap results from a collective effect. 
In the KNM [y], as the coupling is decreased, the RKKY correlations become larger in 
magnitude. The situation is also similar here. Fig. 2a displays the magnetic structure factors. 
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Sf,c{^) — j/^i,m < '3f,c,i^f,c,m > exp[z(/ — m)k], of localized and conduction electrons 
respectively. These structure factors have a maximum at 2ki?c= tt for = 0.2. When 
the coupling is increased, the cusp flattens out, indicating the vanishing of the RKKY 
correlations. We emphasize that even in the RKKY regime, there is no long-range order 
or even a power-law decay of the correlations because the tendency to magnetic order is 
thwarted by the gap induced by the Kondo mechanism. Finally, we note that our results 
are in agreement with those of Igarashi et al. fl^ who also studied the effect of lattice of 



local moments on a Heisenberg chain using exact an diagonalization method. 

Away from half-filling, the situation is less clear because the spin configuration of the 
ground state is not known. Actually, this information is important in a numerical calculation 
since it allows a significant reduction of the Hilbert space under study. There are however 
two limits where one can gain some insight. The first one is the situation with only one 
conduction electron, which is identical to the KLM with one electron. The theorem of 
Sigrist et al. [|Tl|] applies: the ground state is an unsaturated ferromagnet having a total spin 
St = {N — l)/2. The second one is the limit of infinite coupling, in which all the conduction 
electrons form on-site singlets with the localized spins, so that N — Nc spins remain free. The 
ground-state has 2(^-^'=)-fold spin degeneracy. In this limit, the degrees of freedom consist 
in local singlets and unpaired / (up or down) spins. The spin degeneracy can be partially 
lifted by the introduction of the hopping term. The Heisenberg term has no role because 
there are no free conduction electrons. This situation is identical to the strong coupling 



limit of the KLM, where Sigrist et al. have demonstrated that the ground state is a 
ferromagnet with St = {N — for all p 7^ 1. Their result also applies in the present 

case. In order to clarify the behavior of the system in the whole range of parameters, we 
have performed an exact diagonalization for lattice sizes of 4, 5, 6 sites with open boundary 
conditions. We have found that there is a critical value of the Kondo coupling J^? which 
depends on the band filling. When Jk > Jc-, the ground-state is always a ferromagnet at all 
p 7^ 1. Below this transition point, the spin configuration of the ground-state evolves from 
the maximum value of the spin to the minimum. There is a transition region where the spin 



of the ground state can have an intermediate value. But at smaller J^, finite size effects 
become important and trends are less clear. We display in Fig. 3 the energy difference 5E 
between the triplet and the singlet lowest energy states for = 4, 5, 6 sites with Nc = 2, 3, 4 
electrons respectively. For the 4 sites case, 5E \s always positive below the transition point, 
while for 5 and 6 sites its sign is not constant reflecting the flnite size effects. We also present 
the interesting case of 2 electrons in 6 sites. It can help us understand the evolution from 
the maximum spin ground state S't = 2 to the minimum St — 0. At strong coupling, the 
ground state has St — 2, there is a narrow region near the transition point where the spin 
configuration of the ground state is = 1. Then it becomes = at smaller couplings. 
One can also observe that when the density is reduced, the value of Jc decreases. Having 
in mind these results, we use the DMRG to compare the lowest energy states with St — 
and {N — Nc)/2. Since in the DMRG, the states are labelled by only the 2;-component of 
the total spin, at each value of the couphng, we have calculated the lowest energy of the 
states with S^ — 0, ±1 and rk{N — Nc)/2, such that the last number is an integer. When 
these energies are all different, we have found that the lowest value is that of the state with 
Sj: — 0. It is thus reasonable to think that the ground-state is a singlet. This occurs when 
Jk < Jc- But when Jk > Jc all the states have the same energy ( in the order of 10~^) and 
we conclude that the ferromagnetic state is preferred. These findings corroborate the exact 
diagonalization calculations, namely that there is a cross-over from a singlet ground state to 
a ferromagnet at the densities that we have investigated. We caution that in the finite size 
study there is a transition region near where the ground state can have an intermediate 
value of the total spin. A detailed study of this region is left for the future. 

Let us now reexamine the magnetic structure factors of the / electrons shown in Fig. 2. 
In the strong coupling limit, Jk = 2, starting from half-filling (Fig. 2a) where each / electron 
is screened by a conduction electron, the system is in non-magnetic state. The reduction of 
the density is equivalent to the insertion oi N — Nc free spins within Nc singlets. The problem 
is analogous to the insertion of holes in the strong U limit of the Hubbard model. At p = 0.9 
(Fig. 2b), we found that the ground-state has St — {N — Nc)/2, thus, contrary to the strong 
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U limit of the ID Hubbard model, the ground state is an unsaturated ferromagnet; this 
extends the results of exact diagonalization to larger systems. Although the system is in the 
ferromagnetic state, the cusp at /c = is very small, we believe that this is due to the fact 
that the density of the singlets is still much larger than that of the unpaired / electrons, so 
that ferromagnetic correlations are hidden since the effective coupling between the unpaired 
/ electrons is small [|12|. As one decreases the band filling (Fig. 2c, 2(i), the maximum at 



= steeply rises. The lower is the density, the higher is the maximum, since the number 
of the ordered spins is greater. At intermediate couplings, the FM state is destroyed for 
p = 0.9 around Jc = 1.2 and p = 0.7, around Jc = 1. One can clearly see that there is 
no dominant feature in Sf{k) in these two cases when = 1. But the FM state is still 
stable at p = 0.5 as indicated by the maximum at k = 0, the transition occuring only at 
Jc = 0.8. Thus, the FM — PM phase boundary is shifted towards smaller Jk when the 
band filling is decreased. This finding is consistent with the existence of ferromagtism for 
all Jk > of the case with one electron. In the small coupling regime, the ground state 
is always a singlet. An incommensurate peak appears in Sf{k) at 2kpc = Trp at the three 
densities studied, showing that the RKKY correlations are dominant. 

In the weak region we see the ground state is still a singlet despite Nc < N . Nozieres 
|§] has argued that at low temperatures, there will not be enough conduction electrons 
to screen out all the localized spins. This incomplete Kondo screening could lead to an 
ineffective RKKY interaction. We see that in our case the nature of the screening process 
seems to be very different from that of the one-impurity case since a spin-compensated state 
exists even though < N . A scenario discussed in [jl3[ , to avoid the Nozieres exhaustion 



problem, was that a significant part of screening is done by the / electrons themselves. 



through the action of intersite RKKY correlations. Following Blankenbecler et al. the 
condition for the quenching of a spin located at the origin gives rise to the compensation 
sum rule: 

< ^/.o >= - E < shs'i, > - E < si,sio > (2) 

1^0 I 



We have verified this relation for Jk — 0.5 where a Kondo state is expected since the 
ground state is a singlet and the magnetic correlations are short ranged . The left hand side 
of (2) is always equal to 0.25. The first term of the rigth-hand side of (2), which gives the 
contribution of the / electrons to the screening slightly increases when the density is reduced. 
It takes the values 0.08, 0.08, 0.11 and 0.12 for p = 1, 0.9, 0.7 and 0.5 respectively when the 10 
first neighbors are included. At the same time, the contribution of the conduction electrons 
which is measured by the second term decreases. It is equal respectively to 0.12, 0.11, 0.10 
and 0.08. Even at half-filling where the number of conduction electrons is equal to that of 
the local moments , the "self-screening" is non negligible. However, when Jk is increased, 
strong on-site singlets are formed and this screening mechanism is no longer possible because 
RKKY correlations between / electrons are suppressed. The major contribution to the 
compensation sum rule comes from the on-site electron-spin correlation — < S^qSJq > 
which is displayed in Fig. 4. This quantity is found to be very close to 0.25p in the strong 
couphng regime. Therefore, in the large Jk limit (2) is satisfied only at half-fiUing, since 
the / electrons are FM ordered when p I. Thus the depletion effects occur because the 
formation of a coherent singlet ground state necessitates the action of the antiferromagnetic 
RKKY correlations. 

So far, we have mainly discussed the action of the electrons on the lattice of local mo- 
ments. We now wish to consider the inverse problem. For the conduction electrons, we 
have found that Sc{k) is fiat at strong and intermediate couphngs for p — 0.9 (Fig. 26). We 
note that for p = 0.5 and p = 0.7, a small peak arises at /c = 0. This is a consequence 
of the /-spin induced ferromagnetic correlations of the conduction electrons (Fig. 2c, 2d). 
Dominant structures appear in the structure factor Sc{k) at 2kF^ only for small Jk- 

In Summary, we have studied a model Hamiltonian which describes the physics of mag- 
netic moments coupled to strongly correlated electrons in ID. We have found that at half- 
filling, the system is a spin gaped insulator for all non zero values of the Kondo coupling. 
When the system is doped, a FM state is stable at strong coupling. We have argued that it 
is a consequence of the depletion effects. At weak coupling, the ground state is a collective 
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singlet, we have shown that a significant part of the screening is made by the local spins 
themselves through the RKKY interaction. Finally, we have reached the same kind of con- 



clusion in a recent study of the ground-state properties of the one-dimensional KLM . It 
seems that the Hamiltonian (1) and the KLM belong to the same universality class. Indeed, 
in the KLM, there is no explicit interaction between the conduction electrons. However, it 
is well-known that the Kondo exchange term induces an effective interaction between these 
electrons. It is thus likely that the Hamiltonian (1) and the KLM are not fundamentally 
different. 
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FIGURES 

FIG. 1. The singlet triplet gap A (circles) and the binding energy Eb (diamonds) per site 
versus the Kondo coupling. 

FIG. 2. The magnetic structure factor of the localized spins and of the conduction electrons 
(inset) for (a) p = 1, (6) p = 0.9, (c) p = 0.7, (d) p = 0.5 at Jk = 0.2 (circles), Jk = ^ (diamonds) 
and Jk = (stars). 

FIG. 3. Energy difference 6E between the lowest states with S't = 1 (or 2) and S'r = for 
systems of 4, 5 and 6 sites with 2, 3 and 4 (or 2) electrons respectively. In the inset, 5E between 
the lowest states with S^ = Q and = {N - Nc)/2 for p = 1/2. 

FIG. 4. The on-site electron-spin correlation versus the Kondo coupling at p = 1 (circles), 
p = 0.9 (diamonds), p = 0.7 (stars) and p = 0.5 (triangles). 
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